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Nonlinear oscillations within a plasma slab are analyzed. Two types of solutions are 



^ ■ found, depending on the initial value of the electron density. The first represents regular 

^ ■ oscillations within the plasma slab, while the second gives rise to explosive growth at the 

^ \ slab centre or at the edges. The results are discussed. 

ly^ , PACS numbers: 52.35.Fp (Electrostatic waves and oscillations), 52.35.Mw (Nonlinear phenomena: 
U • 

' ^ ' waves, wave propagation, and other interactions) 

^ ■ 

' — I. INTRODUCTION 

> ■ 

O I A few plasma physics problems can be solved exactly (e.g. Akhiezer and Lyubarskii, 1951; 

O . Dawson, 1959; Davidson, 1972; Shivamoggi, 1988; Stenflo and Yu, 1997). The plasma distur- 

O I bances have in those cases been allowed to be so large that no expansion techniques are applicable. 

^ ■ The new exact solutions can be useful in interpreting observed large amplitude wave phenomena 

_ ^ ! in detail, as well as for verifying new approximation schemes and numerical methods in the study 

>~>' of nonlinear effects. We shall in the present paper reconsider one of the simplest possible situa- 

Ph' tions, namely that of one-dimensional (dy = d^ = 0) oscillations in a cold plasma slab {—d <x<d) 

. !_H I where the electrons are mobile, whereas the ions are immobile and form a fixed background with 

X' 

' the constant density hq. The special case where the electron density n is only a function of time 
and where the cold electron fluid velocity v{xj) satisfies the boundary conditions v{±dj) = 
has been solved previously (Aliev and Stenflo, 1994; Stenflo, 1996). Here we shall use the same 
boundary condition but consider a more general and useful class of solutions for n{xj). 
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II. BASIC EQUATIONS 



Our governing equations are 



<9,n + <9,M=0, (1) 
dtV + vd^v = -{e/m)E, (2) 
d^E^{e/£Q){nQ-n), (3) 

where E is the electric field, e is the magnitude of the electron charge, m is the electron mass, and 
Co is the vacuum dielectric constant. By eliminating the electric field in Q and Q we immediately 
find the equation 

d^{dtv + vd^v) = -CofX'^ -n/no), (4) 

where cOp = noe^/eom. 

We shall below study a particular solution where n{xj) and v{xj) satisfy the system of two 
coupled equations (U]) and ©, and where v{±dj) = 0. Our boundary condition means that the 
electrons are always contained within the slab, and that j'^^n{x,t)dx accordingly is a constant 
(= 2d no). By means of straightforward calculations it is then easy to verify that a solution of ^ 
and © is 

n(xj) = r— ; -, (5) 

and 

v{x,t) = -dC0pK{y)sm{(0pt), (6) 
where y{x, t) is given by the implicit relation 

x = d\y+K{y)co%{(Opt)\, (7) 

and the prime denotes differentiation with respect to y. Here K{y) is an arbitrary function satisfying 
K"(±l) = 0. We shall in the present paper choose the simple function K{y) = ay{l —y^), where a 
is a dimensionless constant. With this choice of K we note that dTJ) is a third order polynomial in 
y, and that an analytical expression y{x, t) thus can be presented by solving 

x = d\y + ay{l-y^)cos{(Opt)]. (8) 

However, here we do not give the explicit expression for y{x,t). Other choices of K lead to more 
elaborate expressions for y{x,t) and will therefore not be considered. Corresponding to (El), ^ 
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FIG. 1 : Nonlinear periodic and symmetric density oscillations. To the left the normalized density n is plotted 
and to the right the normalized velocity v is plotted, both as functions of the normalized length xfor different 
normalized times t, with a = —0.4. At t = the maximum density can be seen at x = 0. When a > we 
have the same dynamics shifted by a quarter of a period. 

and Q now reduce to 

"^""'^^^ l + a(l-3/)cos((0/)' 

and 

v{x,t) = -d(Opay{l - y^) sin{Q)pt) . (10) 

We note that n(±J,0) = no/(l — 2a). With the normalization n/nQ,v ^ v/dcOp, x/d, 
t (Opt, and J ^ 1 we plot the density and velocity profiles for different times. One can then verify 
that the total number of electrons is conserved during the oscillations. First, when 0<a<l/2we 
show (in the left panel of Fig. 1) n{x,t) for different normalized times t, using a = —0.4. We note 
the symmetric and oscillating character of n. The corresponding velocity profile v{xj) can be seen 
in the right panel of Fig. 1. For these regular oscillations, the sign of a characterizes the initial 
value of the electron density n{x,0): when a > we start with higher edge density at x = ±d, 
while for a < the density has an initial spatial maximum atx = 0. 

The denominator in Q can approach zero at a finite time to, resulting in explosive growth. Due 
to the symmetric character of the density distribution the explosion for t = to will either be at (a) 
the edges of the slab (x = ±d, i.e. j = ±1) or (b) at the centre of the slab (x = 0, i.e. y = 0). In 
the case (a), the explosion time to is given by cos(cOp?o) = l/2a and |a| > 1/2, while in case (b) 
the explosion time is found from cos{(Opto) = —1/ a and | a| > 1. Thus, in case (b) we can have 
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FIG. 2: Explosive instability and wave breaking. To the left the normalized density n is plotted and to 
the right the normalized velocity v is plotted, both as functions of the normalized length x for different 
normalized times t, with a = —0.8. Att = the density is focused but finite atx = 0. As the system evolves 
we can see wave breaking and density spikes developing at the edges of the slab. Here the case a> Ois no 
longer associated with a simple phase shift, since wave breaking then occurs att = 0. 

explosive growth both at the edges and at the centre, depending on the initial conditions, while this 
it not necessarily true for case (a). The explosive instability occurs in the wave-breaking regime, 
and our cold plasma model thus breaks down at this stage (see Fig. 2). A change of sign of a shifts 
the density distribution by half a period. Note however that a > gives rise to unphysical initial 
conditions, i.e. n(±rf,0) < 0. 



In the present paper we have generalized the work of Aliev and Stenflo (1994) to allow for 
a class of solutions where the electron density n{x,t) is an explicit, and comparatively simple, 
function of x. We have presented figures which show the behavior of the solutions for different 
initial values n{x,0). From those solutions we can very clearly see that n{xj) has an oscillating 
behavior when < |a| < 1/2, while if |a| > 1/2 we have explosive growth, i.e. n grows like 
l/{t — t())^ where the explosion time to is of the order \/a(Op. We think that the present solution 
can be generalized to cold many-component plasmas (Amiranashvili et al., 2002), to plasmas 
where no is a function of x (Stenflo and Gradov, 1990; Karimov, 2002), to cylindrical plasmas 



III. CONCLUSION 
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(Stenflo and Yu, 1997; Karimov, 2005), to Pierce beam plasma systems (Matsumoto, Yokohama, 
and Summers, 1996), as well as to other kind of plasmas (Vladimirov, Yu, and Tsytovich, 1994). 
Alternative solutions are of course also possible (Polyakov, 1995). 

Our governing equations are easily generalized to the relativistic case by replacing Eq. ^ by 
dfp + vdxp = —eE, where p = mv/ a/1 —v^/c^, c is the velocity of light, and m is the electron rest 
mass. At present we have however not been able to generalize Q and Q to the relativistic case, 
and we have therefore limited our presentation to the regime dcOp ~ v <^ c. The opposite limit is 
of course also of much interest, but it has to be delegated to future numerical work. 
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